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Abstract 

Effective thermal masses of bosonic particles in a plasma play an important role in 
many different phenomena. We compute them in general supersymmetric models at 
leading order. The origin of different corrections is explicitly shown for the formulas 
to be applicable when some particles decouple. The correct treatment of Boltzmann 
decoupling in the presence of trilinear couplings and mass mixing is also discussed. 
As a relevant example, we present results for the Minimal Supersymmetric Standard 
Model. 
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1 Introduction 



If supersymmetry is realized in nature it would have manifold and very important impli- 
cations on the history of the early Universe. In fact, much effort has been devoted to the 
study of supersymmetric cosmology and truly supersymmetric solutions to old cosmological 
problems have been proposed (while new problems have also arisen. See [Q for review and 
references). However the (weak scaleQ) supersymmetric generalization of the Standard Model 
(SM) is not uniquely defined. First of all, the introduction of arbitrary soft supersymmetry 
breaking parameters, to prevent the mass degeneracy between ordinary and supersymmetric 
particles, generates a lot of freedom and second of all there are various options related to 
the particle content and the gauge group definition. This limits the generality of the pre- 
dictions that can be made although still permits to confront different classes of models and 
theoretical assumptions by examining their cosmological implications. 

At the high temperatures of the early Universe, supersymmetric particles would be ther- 
mally pair created and would populate the plasma. One of the simplest consequences of this 
fact is that the effective thermal mass of a generic particle immersed in that plasma would 
be changed due to interactions with supersymmetric ambient particles. It is obvious that 
knowledge of these effective thermal masses is fundamental to describe the behaviour and 
properties of the plasma. Moreover, it is well known that these quantities play a crucial role 
in many interesting aspects of the evolution of the early Universe. Various examples follow. 

In the case of gauge vector bosons (see e.g. [fj, U) the effective thermal mass for longi- 
tudinal components corresponds to the usual Debye mass, i.e. the inverse screening length 
of electric potentials in the plasma. At leading order (one loop in perturbation theory) it is 
ttid ~ gT, where g is the corresponding gauge coupling constant. Transverse components 
have instead zero thermal mass at leading order. For abelian gauge bosons this is true also to 
all orders, corresponding to the non screening of magnetic fields, but for non-abelian gauge 
bosons a magnetic mass of order g 2 T is expected to appear non-perturbatively. Supersym- 
metric particles in the plasma will have an influence on Debye masses (see e.g. |4j]). In this 
paper we will consider only thermal masses at leading order so that magnetic masses will be 
taken to be zero. 

Let us turn now to scalar thermal masses. As pointed out by Kirzhnits and Linde ||, 
spontaneously broken symmetries are generally restored at high temperatures (see also 0). 
This can be understood in terms of the effective thermal mass of the (Higgs) scalars driving 
the symmetry breaking. Consider as a particularly relevant example the electroweak gauge 
symmetry. Call <fi the Higgs field responsible of the breaking. The one-loop approximation 
for the effective potential of 0, including the effects of finite temperature is of the form 

T) = ^T 2 - m 2 )0 2 - ET^f' 2 + ^A(T)0 4 , (1) 

where E, k and A(T) are some functions of the masses and couplings, easily calculable in a 

1 We concentrate here on temperatures of that order, relevant for example in studies of the electroweak 
phase transition. 
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given model. For low temperatures the negative T = mass squared dominates favoring the 
formation of a condensate, while at sufficiently high T, the (leading order) Higgs thermal 
mass a/kT ~ [coupling] x T dominates over the negative T = mass disfavoring a non-zero 
condensate. Furthermore, it is clear that knowledge of the thermal mass for the Higgses 
allows an estimate of the critical temperature of the transition (T c 2 ~ m 2 / k). 

The order of such transitions is also related to the value of thermal masses in a more 
indirect way. From ([I]) it is clear that the presence of the non-analytic cubic term causes the 
transition to be first order. In fact, the jump in the order parameter is 

T c X(T c y ^ 

[here T c is defined by the coexistence of two-degenerate vacua in ([[])]. The quantity (H) 
is of the utmost importance for the viability of electroweak baryogenesis (for review and 
references see e.g. [[?)). Now, the cubic term in ([!]) is a purely finite temperature effect and 
comes from the interaction of the Higgs field with the static modes of different species of 
bosons in the plasma (fermions do not contribute to this term because they do not have 
static modes). In fact, each bosonic degree of freedom, with (T = 0) field dependent mass 
Mj(0) contributes to the potential V(<p,T) a term 

A t V = --^iM?(<p)f 2 . (3) 

Beyond the one- loop approximation for the potential, every mass Mj(0) in ([3]) should be 
substituted by the corresponding effective thermal mass, obtained from 

M^^M^ + ^-T 2 , (4) 

where the last piece comes from the interaction of the particles with the surrounding plasma. 
Substitution of (|4]) in @ resums an infinite series of higher order diagrams, the so-called 
Daisies. The net effect of this resummation is to screen the cubic term in ([!]), effectively 
reducing the E parameter and thus weakening the strength of the phase transition. In the 
Standard Model, where the dominant contribution to the cubic term in the potential comes 
from gauge bosons, the screening of the longitudinal modes is very effective while it is zero at 
leading order for the transverse modes. Then, daisy improvement of the effective potential 
leads to a reduction of the strength of the transition || roughly by a factor 2/3. In the 
Minimal Supersymmetric Standard Model (MSSM) if stops are light they give the dominant 
contribution to the cubic term in the Higgs potential (see || for the effect of Debye screening 
on the electroweak phase transition in the MSSM) and the final strength of the transition 
will be sensitive to the value of stop thermal masses. 

Besides the effects explained, (bosonic) thermal mass corrections are very important 
because they represent the starting point of a resummation of perturbation theory (see e.g. 

|HJ). Such resummation is necessary to take care of the infrared problems that plague 
theories at finite temperature if they contain massless bosons in the symmetric phase, e.g. 
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Yang-Mills theories [ XT | . The problem appears when we probe our system to low scales^ 0{gT) 
compared with the temperature T. At this scale, an infinite number of diagrams can give 
contributions of the same order and to improve the usual perturbative series they need to be 
resummed. The effective thermal masses provide then an IR cut-off taming the perturbative 
expansion^. One example is provided by the cubic term in the potential discussed previously 
Its non-analytic behaviour signals its infrared singular origin: it comes from (bosonic) zero 
Matsubara frequency modes. Note that fermions do not cause infrared problems because 
they do not have zero Matsubara modes. In fact, at sufficiently high temperatures (or for 
distances much larger than 1/T), fermions decouple from the effective 3D theory at finite T. 
For that reason we concentrate here on bosons only. 

Other examples where effective thermal masses play a role (in supersymmetric contexts) 



are: studies on the 3D reduced effective theory in the MSSM [[12], analysis of charge and 
color breaking minima [I3| at finite temperature |TJ]], non-restoration of symmetries at very 



high temperature in general supersymmetric models ||15|| , inverse symmetry breaking at 
some range of temperatures [M-[|T7|], different details of the spontaneous mechanism for 



electroweak baryogenesis [T^j, etc. 

The aim of this paper is then to compute thermal masses for bosons (scalars or gauge 
vectors) in general softly-broken supersymmetric models (section 2). In subsection 2.1 these 
masses are presented for temperatures much larger than all particle masses. In that case all 
particles in the theory are thermally produced and form part of the plasma. In subsection 
2.2 we present the more complicated case in which the temperature is lower than the mass 
of some particles which decouple from the thermal plasma and then do not contribute to 
the effective masses of other particles. Section 3 applies these results to the particular 
case of the MSSM (some of the results presented have already appeared in the literature 



ni 13 P). 



2 General Softly-Broken Supersymmetric Model 



Since Bose-Einstein and Fermi-Dirac distributions are different, the thermal bath is pop- 
ulated by different amounts of on shell fermions and bosons. In such a sense, in a SUSY 
theory, temperature effects can invalidate various cancellations implied by the symmetry 



between fermions and bosons f22|| . This observation is very relevant in particular for the 
computation of effective thermal masses. 

As is well known, only self-energy diagrams which are quadratically divergent at T = 
contribute to the leading thermal masses. Typical diagrams that enter such calculation are 
depicted in figure 1. Although the second diagram is not quadratically divergent it can give 
a contribution in the presence of Boltzmann decoupling and should be kept. Note that for 



2 Here g stands for a typical gauge coupling or a Yukawa coupling. For power counting quartic scalar 
couplings are A ~ g 2 . 

3 Of course, there remains an infrared problem for transverse gauge bosons, associated with physics at the 
scale g 2 T. 
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our purpose external momentum can be set to zero. If fermion-boson cancellations were still 




Figure 1: Different types of diagrams contributing to thermal masses and responsible for the indicated 9 
symbols. 

operative at non zero temperature in the supersymmetric case we would obtain zero thermal 
masses. However, it can be shown that fermionic contributions come with an extra factor 
(—1/2). More explicitly, if a bosonic integral gives 

/ ft = /<A 2 + T 2 ) + ... (5) 

where A 2 is the T = quadratic divergence and kT 2 the associated finite temperature 
contribution to the thermal mass, the fermionic counterpart will be 

If = -^ 2 - l -T 2 ) + ... (6) 

Then, instead of cancellation of thermal masses there is a reinforcement: 

I b + I f = ^T 2 + ... (7) 

Explicit examples of this effect can be found in the next sections. 



2.1 Thermal masses in the limit T ^> M 



As we will see, in the limit that the temperature is much larger than any mass in the theory 
the contributions to the various self energies depend only on the gauge structure of the 
theory and the dimensionless parameters of the superpotential W which reads 

W = + yWijk^jfa. (8) 

Latin indices i,j,k,... will be used for scalar fields. The corresponding fermionic partners 
carry a tilde: k,l,.... Latin indices a, b, c, ... are reserved for gauge bosons and the tilded 
version for gauginos. Unless stated otherwise sum over repeated indices is always implied. 
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The full scalar potential is then: 



E 



dW 



1 2 
+ ~ 2 9l 



+ m 2 |0i| 2 + -Bijfafa + —Aij k (j)i(j)j(j) k . 



The remaining soft breaking termsQ are gaugino masses: 

KM,K + h.c. 



(9) 



(10) 



Leading order thermal masses for scalars can be obtained simply taking derivatives from 
the one-loop finite T effective potential which reads 



v{<p) = v {<t>) + v l {<p) + v T {4>), 

where Vi((f>) is the T = one- loop correction and 



'.ID 



v T (<j>) 



rp4 



2 T,9sJs 



J"2 



with J s = J + (J_) if the s th particle is a boson (fermion) with g s degrees of freedom (defined 
negative for fermions) and 



POO 

J±(y 2 ) — I dxx 2 \og 
Jo 



(±)e 



-y/x 2 +y 



(13) 



The behaviour of J±(m 2 /T 2 ) is very simple in two different limit cases. First of all, 
the expansion of J±(m 2 /T 2 ) for large values of m/T give contributions that are exponen- 
tially suppressed ~ e~ m//T , while the expansion for small values of m/T gives the leading 
contributions 0(T 2 ): 

V t-^H 9 b m 2 + — 9f m ) + - (14) 

boson fermion 

Then, we will simply use a step approximation for the effective potential to compute the 
thermal mass corrections: 



rp2 rp2 

V t = 24 S 9bml6 b + — 9fm 2 f 6 f 



boson 



18 



(15) 



fermion 



where the sum is on all mass eigenstates calculated in the theory at zero temperature and 
9bj = 1 if rribj T and if rribj ^> T. Of course this is a crude approximation but gives 
the correct results in the two limiting cases of interest. Now we study the limit in which the 



4 If the model contains matter fermions in the adjoint representation of some group, soft masses that mix 
them with the corresponding gauginos can be written. However these soft terms are generically absent in 
super gravity scenarios. 
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temperature is the larger mass scale. In such a case 9bj = 1 for all bosons and fermions. The 
analysis is very simplified due to the fact that we can forget the complications induced by soft 
breaking mass terms and supersymmetric massive parameters present in the superpotential. 
Setting then all #'s to 1 in ( Jl5| ) and using the fact that 

StrM 2 ((f)) = 3TrMy + TrM 2 s - 2TrM 2 F = £ g s m 2 s ((j)) = K- 2g a D a TrT a } (16) 



with K a field independent constant and D a = (f)*T^(pj, we obtain 



Vr 



T 
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m, 



D — term + K' . 



(17) 



The constant term is irrelevant for our purposes and the D-term vanish if we further assume 
TrY = for any U{1) gauge group present. 

Then, up to some irrelevant constant, we can write 



V T (4>i,T) 
And from this, 



12 

16 



12 9 f 



r- 



i,k 



i.k 



d 2 V T (<P,T) 
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Y.w lkl w* kl + 4j2g 2 a (T a T a 



1 ■>■ j 



l k,l 



(19) 



Writing {T a T a ) i j = C a (R)5ij and using a convenient basis for the fields 0j we get 



T 2 



k.l 



(20) 



which gives the thermal mass corrections for scalars. This diagonal correction should be 
added to the T = mass matrix. The eigenvalues of this thermally corrected matrix are the 
end point of our calculation. 

The leading order thermal masses for longitudinal gauge bosons, Hy, get contributions 
from scalar, fermion and gauge boson loops plus their supersymmetric partners (note that 
we can describe the chiral supermultiplet contributions either as S + S or F + F . We use 
both below in the understanding that no F corresponds to any S) 



n 



, = n<? + nf + uP + uP + up + uP = 6 -uP + 3 n^ + fn^, 



(21) 



where the last equality follows from supersymmetry as explained above. Now, the vector 
contribution is 



17(1) : UP = 0, 



SU(N) : UP = jg 2 N T 2 



(22) 



6 



and according to our rule gaugino loops contribute half this result. The contributions from 
scalars and (chiral) fermions are 



6 S F 



(23) 



with Tr(T a T h ) = t2(R)5 ab . Sfermion loops contribute twice as much, while higgsinos give 
only a half. The contribution from non-chiral fermions is twice larger than that from chiral 
ones. 

The final result is then 



n 



U(l) 



n 



SU(N) 



E^ 2 + E^ 



En 

A 



N + ^2h(Rs)+J2t 2 (R F ) 

S F 



;9nT 2 



N + ^t 2 (R A ) 

A 



(24) 
(25) 



where the index A runs over chiral supermultiplets. 

We see explicitly that all possible self energies depend only on the gauge quantum num- 
bers of the spectrum and on the Yukawa couplings that appear in the superpotential. 
In practice, at very high temperature the masses of the underlying T = are irrelevant. In 
such a case, the computation of the leading thermal corrections is simplified and they can 
be derived directly from an exactly conformal supersymmetric theory. 



2.2 Thermal Masses for general T 

The study of the case in which the scale of the temperature is not the dominant one is 
a bit more involved. The mass scales present in the theory, aside from possible non zero 
background fields, are the soft susy breaking terms and the massive coefficients in the bilinear 
terms of the superpotential. These scales can have very different values and there always 
exist some range of temperatures in which decoupling and mixing effects have to be taken 
into account. 

The obvious new effect is the Boltzmann decoupling of particles of mass m ^> T and 
then of their contribution to thermal masses of other particles. This effect will be taken into 
account by writing every contribution with the corresponding 6(T — m) that will take care 
of the decoupling in a step approximation. 

As long as field-background effects can be neglected, i.e. as long as no particle is de- 
coupled because of a large background dependent mass, the thermal self-energies will only 
mix particles with the same quantum numbers. The reason for this is that leading thermal 
masses arise from quadratically divergent diagrams that can already be drawn in the T = 
unbroken theory. As an example consider W3 — B mixing. One can certainly draw diagrams 
at T = that mix those particles at one-loop. However when summing such diagrams over 
complete SU(2) multiplets these contributions cancel. At finite T this is reflected in the 
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following non-diagonal self-energy: 

U WslBl = ^ gi g2T 2 [Trs(9YT 3 ) + 2Tr F (9YT 3 )}, (26) 

which gives zero when all #'s are 1. In a background that breaks £77(2) x £7(1) if some 
particle acquires a mass larger than T its contribution to (|26| ) will drop and only then a non- 
zero contribution will result. Here we will always assume that background masses are always 
smaller than the temperature (which is usually the case in most applications of interest) 
so that we will not encounter this complication. In that case one can compute thermal 
masses at zero background (corrections from non-zero background effects will be suppressed 
by powers of T). The prescription to obtain thermal corrected masses is then to write the 
T = mass matrices in whatever field background one is interested (provided it is smaller 
than T), add the thermal corrections and afterwards rotate or diagonalize the mass matrix. 

Setting then zero background we can in principle compute thermal self energies using an 
interaction basis or a mass eigenstate basis. The first option is more convenient and it is 
simple to rotate to the mass basis in particular cases (note that 0's are naturally defined in 
the mass basis, so that, to decouple some particle the rotation should be made). We will 
express our general results in terms of some convenient 6 symbols which vary with the origin 
of the contributions as shown in figure 1. The rules to rotate these symbols from one basis 
to another are explained below. 

The fields in interaction basis, fa, can be written as a linear combination of the mass 
eigenstates <p a : 

fa = U?<p a , (27) 

where we stress that the unitary matrices U diagonalize the M 2 mass matrix calculated at 
zero background and zero temperature. The symbol 9y comes from the contraction of <p* — (f)j 
to close the loop as shown in figure 1. It is defined by rotating to the mass basis as 

% = Qu? Va I7? VB = U i* d °0 U j = U i* 0aa U h ( 28 ) 

with 

9 a (3 = 8aa3a(3 = I n • j " _ rp (29) 



The 6{j symbol defined applies both to fermion or boson contractions. 

We define also the 4-index symbol 9fj for the second and third diagrams shown in fig. 1. 
For these objects the rotation from the interaction basis to the mass basis is 

e» = u? uf 0% ur vf. (so) 

But now there is a difference for the fermionic and bosonic case. For fermions we have simply 
(tildes omitted) 

QfiP = @aa9f3p, (31) 

while for bosons: 

acta _ @ aa ~ QfiP (nn\ 
U W ~ m 2 _ m 2 ■ 
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The reason for this is the following: note that the bosonic diagram is not quadratically 
divergent (in particular 9gg = if 6 aa = 6pp = 1). However it contributes to the thermal 
masses if one of the particles running in the loop, say a, decouples. In such case, the diagram 
behaves effectively as the first one, with the heavy line in the loop collapsed to a point. In 
other words, in the effective theory that results after integrating out the heavy particle there 
are new quartic couplings proportional to l/m 2 a . The symbol ( |3"2"| ) takes this into account. 

There is another effect we have to mention before presenting the results. Suppose that the 
scalar fields fa and 4>j have the same quantum numbers but opposite abelian charges. The 
mixing fa — <p* by thermal mass effects is not possible in the non-decoupling case analyzed in 
the previous subsection [see eqs. (18) and (19)] but becomes possible in the case that thermal 



contributions from some particles are Boltzmann suppressed. We allow for such possibility 
in our general formulas. The corresponding thermal self-energy will be denoted by IX^*. 

Also note that we give our results in terms of thermal polarizations and 0's for complex 
scalar fields. This assumes that real and imaginary components behave in the same way, e.g. 
they decouple together when some mass parameter is made heavy, etc. This is no longer the 
case in the presence of large backgrounds, which we assume not to be the case, or for singlet 
fields. In this last case, real and imaginary components can have different masses and should 
be treated separately. Our formalism can be trivially generalized to take this possibility into 
account using relations like 

(or i aq\ J ®SS — \ \§S T S T + ^S^ 1 



y2 { fss* — 2 ps r s r — Vs^s*] 



and so on. 

The general results are the following: 
A. Scalars 

A.l Yukawa contribution from fermion loops: 



T 2 

1 QSZ 



n*,* = ^W M W* rs 6%. (34) 



24 

A. 2 Yukawa contributions from scalar loops: 

n*,* = ^w rik w;^ n^* = ^w riJ w; kl e kl ,. (35) 

A. 3 Trilinear contributions: 

Note that these terms are proportional to Q\\ and thus give zero in the limit T ^> M. 
T 2 

^-4>u4>i = 7^7 {AikiAj rs 9fg + 2W iM ^ m W* rs fj, rn (6^ n + 6^ gt ) + W^ kl ^ mi W rsn fi* rj 6 s ^ l (36) 



+ 



24 

2W ik i^ km W nrs fi* r 6^ + 2A m W* ra n m rQ s **i + A m W rsm ii) T et\ + (h.c.,i <-> j)] } 



T 2 

n^,^ = ttt {A m Aj ra 9i*s + 2Wikifii m Wj r sti* n (9 l ^ s + +W* kl ^ ir W^ ns fij m 8^ (37) 
+ [2AkiW jmnf i* mr 9ff + A lH w^ nrHm el n + 2w: kl ii ir w jmn ^* ms ef l + (i <-> j)] } . 
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A. 4 Gauge contributions from fermion loops: 

We change momentarily our notation from fa to N a where N refers to a given rep. of 
the group and a to the group index (not to confuse with mass eigenstate indices). Of course 
no change is needed for C/(l)'s. 



rp2 rp2 



(38) 



The second contribution can be non-zero only if the model contains matter fermions in 
the adjoint representation as discussed in footnote 4. 
A. 5 Gauge contribution from scalar loops: 
The general result is: 



T 2 

u N a ,p p = —g z 



T 2 



ri/v p* = — q 1 



rpa rpa n 
i 7Q J <5/3 t/ '5*7 



For SU (N) with all non-singlet fields in the fundamental rep. 



n 



N a ,Pg 



2A 9 



28 NP T$ a Tr s (T a 9) + 5 aP 9 N ^ - -9 aP 



T 2 



^N a ,p* - -^g 1 



For U(1) Y : 



H*,* = ^ [5 tJ Y/Tr s (Y9) + Y^} ; = ^9^9 

A. 6 Gauge contribution from gauge boson loops: 
The general result is: 

T 2 

Rp a ,Pf3 = -T-gAgBTfaT^Oab, 
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For SU(N), when 9 ab = 5 ab 9 aa : 



n 



For 17(1): 



P a ,P« 



IL 



g 2 N C N {Rp)8 a p9 G - 



T 2 



glY 2 5 i3 9 B . 



(39) 
(40) 

(41) 
(42) 

(43) 

(44) 

(45) 

(46) 



B. Gauge Bosons 

As already mentioned, only longitudinal gauge bosons get a non-zero thermal mass at 
leading order. The following thermal polarizations should then be understood as polariza- 
tions for the temporal components n 00 of the gauge fields Vfi, Vq. 
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B.l Scalar contribution: 



T 2 r 

n a fe = —9a9b [{T a ,T }p 1 6p fjPri + 



T f3a T S~f {pM g N^N 5 M a - @Mfj N$ @N~, M* ^ + h.C. 



(47) 



where in principle two different groups, with coupling constants g^, qb are considered. 
For SU (N) and fields in the fundamental rep. M, P, etc: 
The general result is: 



0M a P a 0p Ci M Ci - MR PX^MXP fj 



a. J "a 



where 9 M = (1/N) J2 a ®M a M a - 
For U(l): 



T 2 



K BL = —glY.YiYi 

ij 



B.2 Contributions from matter fermion loops: 



T 2 



n a b - -^-9A9BT^ /3 T^ s ep i3 f ;ri 9p a f !g . 



B.3 Contributions from gaugino loops: 

rp2 



n 



«/. '~9A9Bfdeafghb^i 



(48) 



(49) 



(50) 



(51) 



where f^ bc are the structure constants of the group A ([T a ,T b ] = if a b c T c ). 



3 Minimal Supersymmetric Standard Model 

In this section we apply our general results to a particularly relevant example, the MSSM. It 
is the simplest supersymmetric extension of the SM and is described by the superpotential 

W = f xH 1 -H 2 + h Uz Q t ■ H 2 Ui + h Di E x ■ QiDi + h Ei E x ■ UE % (52) 

embedded into the SU (3) x SU (2) x U (1) gauge group. We will retain full freedom in the soft 
susy breaking terms in order to be as general as possible. However we will assume negligible 
intergenerational mixing. In such case the only fields that mix at zero background are Hi 
and H 2 . 

The scalar-fermionic soft lagrangian reads: 

Lsoft = A Ei Hi ■ LiEi + A Di Hi ■ QiDi + A Vi Qi ■ H 2 Ui 

m\H x ■ H 2 + [£gMg + h.c] + $>?|&| 2 . (53) 

9 i 
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3.1 Limit of very large T 



The bosonic self energies in the case in which the temperature is the larger mass scale, i.e. 
T 3> /i, A^, m^, m 3 , M are obtained from subsection 2.1 directly as 



Urr 


= n fi = -G 2 T 2 + -o 9 2 T 2 + — q 2 T 2 + + /i 2 n )T 2 
3 3 8 72 4 Di' ' 




- ^ 2 2 , 2 2 2 , 1 . 2 T 2 


n - 


- - a 2 T 2 + — a 2 T 2 + -hA IT 2 


TT 


TT ^ 2^2 , 1 ^,2^2 


TT- 


- 2 9i+ , 




= = ^ 2 T 2 + i, 2 T 2 + ^ 2 T 2 , 




= n ff2± = ^ 2 T 2 + l -glT 2 + ^ 2 T 2 , 


n 5i 








n 5i 


= T^ r • 



(54) 

For if x and H 2 we only keep third generation Yukawa couplings. Also, note that particles 
in the same gauge multiplet receive the same thermal mass correction. 

3.2 Explicit formulas in the general case 

If there is not a defined hierarchy between the scales T, /i, A^, rrii, m^, M, we must apply 
the formulas of subsection 2.2 that obviously have as asymptotic limit , for high T, the 
equations presented in the previous subsection. 

In the formulas that follow we write most of the self energies in the interaction basis. 
Also, we use &i = Qa, Q\ = 9 J J , etc, to simplify the notation. For all fields besides Hi and 
H2 the 8ij functions in the gauge basis are diagonal and coincide with the definition in the 
mass eigenstate basis (at zero background). The treatment of ifi i2 is as follows: as is well 
known, there are three mass parameters in the tree-level Higgs potential of the MSSM: 

V = m\\Hi\ 2 + m\\H 2 \ 2 + m\(Hi ■ H 2 + h.c.) + quartic terms. (55) 

Two of these mass parameters (1711,1712) can be traded by the T = vacuum expectation 
values i>i and v 2 [with v \ + v 2 = v 2 = (174 GeV 2 ) and tan/3 = v 2 /vi\ leaving only one 
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free mass parameter, conventionally taken to be the mass of the pseudoscalar Higgs, m^. 
Then we have two scales in the Higgs sector, v (or Mz) and m^, and the only non-trivial 
case at finite temperature corresponds to Mz <T< uia- When rriA 3> Mz one (linear 
combination) of the two Higgs doublets is heavy (~ tua) and one light (~ Mz)- 

In order to obtain the mass eigenstates at zero background, we can work with the full 
doublets. Diagonalization of the 2x2 matrix defines the mixing angle (3q. In the only non- 
trivial case with uia ^> M z , it's straightforward to see that j3 — > (3 so that, in this limit we 
can define the doublets H (light) and $ (heavy) by the rotation 

~H\ = H cos (3 - $ sin (3, 

H 2 = i7sin/3 + $cos/?, (56) 

where Hi = (— Hi, H®*) T . The doublets H and <3> are the mass eigenstates so that (|56| ) 
is our equation 4>i = U£*(p a in this case. From (j5fj) and using rules (|28|) and ( |30|) we can 
express all 9 symbols for Higgs bosons in terms of Qu and 6$, or equivalently 9{T — Mz) and 
9{T-m A ). 

From ( p5| ) it follows 

9 H o H o cos 2 (3 + 0$o$o sin 2 /3, 
(9h°h° — 9$o$o) cos /3 sin /3, 

— (9 h +h+ — 9$+$+) cos/?sin/3, (57) 

and so on. The rest of 9 symbols are trivial to handle. For squarks remember that gauge 
invariance requires equal soft mass ttiq for Ul and Dl so that 9fy L = 9q l = ^q. 

Also note that although ^'s for gauge bosons will always take the value 1 (because they 
are massless at zero background) write them explicitly. 

We also use 

6Tr s (9Y) = -3^(0^ + 9 Slj - 29~ eR .) + 3(9 H ± + 9 H o) 
j 

- 3(9 H ± + 9 H o) + N c J2(9 0lj + 0& Lj ~ + 2^), (58) 



2Tr s (9T 3 ) = E(*%- ~ ^) + N cJ2(% ~ H h ) + - 9 H ± + 9 H ± - 9 H o. (59) 

Note however, that this last trace would be non-zero only in a 577(2) l breaking background. 

In the formulas that follow the reader can easily check sector by sector that fermionic 
contributions are always half of the corresponding bosonic ones. 

A. SQUARKS 

Thermal self energies are diagonal in color space unless a color-breaking background that 
decouples some contribution is present. As we assume this is not the case the color index 
structure is trivial and is suppressed. 
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n 



6 ys AN r 



+ 6fj + 29g9 UL 



+ -^g 2 2 T 2 [M w ± + 39 W3 + 9 dLi + 29^ + 2Tr s (6T 3 ) + 29^9^ + A9 Bl ^ ± 
1 



+ 



432 



g\T 2 [39 B + 6 + 6Tr s (9Y) + 29 B 9 ULi 



U Ri D R . 



(60) 



n 



1 2 2 iV 2 -l 

" 6 93 ~JnT 
1 



36> 3 + 9 t)L + 2%# Dl 



+ -r^T 2 Q9 W± + 39 W3 + 9 bh , + 26 6l . - 2Tr s (9T 3 ) + 29 Dl 9^ + A9 Ul 9^ ± 



+ j^gfT 2 [39 B + e 5 + 6Tr s (9Y) + 20 5 D J + A ^ + A^ 



432 ' 
12 ^ 



+ — T 2 
12 A 



0H° + 0.D fl . + Qh°0d r . 



(61) 



n 



1 o^oiV 2 - 1 r 
-a 2 T £ 

6 y3 AN. 



+ —g 2 T 2 [129 B + A9 0ri - 6Tr s (9Y) + W § 9 Ur ] + A ^ + A% £ 
1 



+ 12" u 



h 2 T 2 



(62) 



where 



n 



D, 



1 2rr2 iV 2 - 1 
6^ 3 ~1aT 



3 ^<? + + 2 ^9^D R . 



1 



+ 



/l 2 n.T 2 
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D 



T 2 



12 



-h 



u 



\Au P ?e% + \u,\ 2 9% + (A Up fj, + A* Vp ^)9%^ 

\ a d p \ 2 9hc + \h\ 2 9hc - (A Dp fi + A* Dp ii*)9 hc H c* 
Rotating from H 1: H 2 to H, $ as explained, the previous A's can be written as 

9p — 9 hc ~ _ . o 9p — 9<&c 



T 2 

A ir = T2 h ° 



T 2 

A C P = —h 2 P 
p 12 p 



\Ap\ 2 ZZ2 ~ZJ + \ Ap\ 2 ~'l ,,,2 



771 



(63) 



(64) 
(65) 

(66) 
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Here P = U Lt ,D Lt , U Ri ,D Ri ; c = 0, ± and 

Ay = A^ sin (3 + pi* cos (3, 
Ap. — A Di cos (3 + pi* sin f3, 

B. SLEPTONS 



Ajj. = A Ui cos P — pi* sin /3, 
= A Di sin (3 — pi* cos (3. 



n 



i 



TEdoT Q6 W± + 36 Ws + + IB- + 2Tr s (9T 3 ) + 20 v 6^ + 4fl e ^ 



48 

+ -^- A g\T 2 99 B + 3^ - 6Tr 5 (^F) + 60^ + A l 



144' 

+ -h 2 E T 2 
12 El 



(67) 



48* 
1 



= -9 2 2 T 2 [Qe w± + 3e Ws + e^+2e &Li -2Tr s (eT 3 ) + 2e eL e^ + 4e UL e w± 



(68) 



1 

12 



n 



+ A% + A% 



(69) 



(70) 



= ±g 2 T 2 [1M B + 60^ + 6Tr 5 (#y) + 120^ 

where the A's follow the same notation used for squarks and now 

A E . = AE i cos (3 + fi* sin f3, A E . = Ae { sin (3 — pi* cos (3. 
C. HIGGS BOSONS 

n H o = ^ 2 2 T 2 [6^± + 3^a + H o + 20 H ± + 2Tr s (6T 3 ) + 20^0^ + 40 # ±0^ ± 

+ ^fT 2 [99 B + 3^o - QTr s (9Y) + 60^o] + A, (71) 

+ E [N c h 2 Di {e &Li + e &Ri + e DL e DR ) + h% (e gLi + 9 Er% + 0^); 



n ff ± = ^glT 2 [m w± +39 W3 + 9 Ht + 29 H o-2Tr s {9T,) + 29^9^ + A9 m 9^ ± 
+ ^g\T 2 [99 B + 39 H ± - QTr s (9Y) + Q9 B 9^} + A, 



(72) 
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n„ 2± 



-g 2 2 T 2 [69 w± + 39 Wa + 9 H ± + 29 H o + 2Tr s (9T 3 ) + 26 & ±9^ + 40^ ± 



+ -^g\T 2 [99 B + + 6Tr s (9Y) + 60^ 



+ A 2 



(73) 



48" 

-^g\T 2 [99 B + 30 H o + 6Tr s (9Y) + 60^o 



g 2 T 2 Q9 W± + 39 w , + 9 H o + 29 H ± - 2Tr s (9T 3 ) + 29m9^ + 40 6 ±9# ± (74) 



+ A 2 



n 



H°H*> 



_}_rp2 

"48 



(si + 9 2 y)0h»h»* ~ 2 9lO H ± H T> 



A 



12- 



(75) 



with 



n 



th2 

~48 



(^2 + gl)8 H £Br ~ 2 9l e H° 2 H°* 



A 



12- 



A 2 
Ai 2 



E \Nc W9^+^\A D f9^ 



+ h 2 E .\A E .\ 2 6% 



^2 

12 



|2^Li 
Er- 



4e^ c 



+ hl.AE-0% 

a i 1 Ef 



(76) 



As an example of how to rotate 0's and IPs consider the case in which only one (com- 
bination) of the Higgs doublets is light compared with the temperature while the other is 
heavy and Boltzmann suppressed (this limit is realized for a large pseudoscalar mass and 
has been considered at finite temperature in studies of the electroweak phase transition). 
We will concentrate in the Higgs loop contribution to Higgs thermal self-energies only. The 
rest of the terms are trivial to handle. In terms of $#,0$, the off-diagonal thermal mixing 
between H 1 and H 2 [eqs. (|75D,([76|)1 has the form 



n 



48' 

It* 

48 



-T 2 {g 2 2 + 0?)(0ho - M + 2g\(9 u± - $± ) sin/? cos /3, 



(9 2 2 + g\){0 H ± - 0*±) + 2g 2 (9 H o - $ o) sin/?cos/3. 



(77) 



In the neutral sector then, setting 0$o = 0, it is easy to obtain 
n^o = n^o cos 2 (3 + n H o sin 2 (3 + 2n H o^o, cos (3 sin (3 = — (g{ + g\) (29 H o + 9 H ± ) cos 2 2/3. (78 

1 2 2 1 /I Sv 
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It can be checked that this is the correct result noting that the Standard Model result is 



n 



scalar 
H 



-XT 



(79) 



(with the quartic Higgs coupling in the potential normalized to V = |A|if | 4 ) while in the 
MSSM, the quartic self coupling of H, defined by Q56|), is A = \(g\ + gl) cos 2 2/3. 

D. GAUGE BOSONS 

D.l SU(3) C 

1 



^N c e g + 2Y J {e ijL +e bL +e ijR +e D[ 



+ E( 



2N c 8 g 



To simplify the contribution coming from squark loops we have used 

9q + ®\ = ^Qq- 



JO) 



Note however that if we were to rotate the squark basis the expression on the left hand side 
should be used. 
D.2 SU(2) L 



— glT 2 



IW W ± - 29 gh + 89 w± + N c ^2(29 ^ + 20 b + 9 Ul . + 9 Dl . ) 
+ ^2(29 PL .+29e L .+9 VL .+9 eL .)+2(9 H o + 9 H ± + 9^ + 9^) (82) 



+ 9 m + 9 6 ±+9 s ±+9 6 



2 n 2 



jjglT 2 [3(9 Ws + 9 W ±) + 129 Wa 9 w ± - 29 gh + S8^ ± 8 



24 



h L . + 6 d l .) 2 + 26 u l Qd Li 



w 3 



+ (0 H o + 9 H ±) 2 + (9$o + 9$±) 2 + 26 6 o8 s ± + 29 k ±9 k0 



(83) 



Here 9 g h gives the ghost piece and we have already rotated the Higgs contributions to 
the H, $ basis using: 



J4) 



and similarly for the charged #'s. We have simplify further our expression using a relation 
similar to flSlD. 
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D.3 U(1) Y 
1 



18 + % + 40 iaj l+SE^t #e tj + 4fl e „ . ) 

i i 
+ 18(0 H ± + fl H o + 0$± + 9 0O ) + 9{6 A ± + 6 & o + + 6 & o) (85) 

+ 2iV c 5:(% L + 9 bL + 16^ + 40 Sa .) + N c Y,{9v + DLj + m v + 46 D 

3 3 3 3 . J J ■> ■> 

j j 

Here, contributions from scalars, and in particular Higgs bosons, have been treated in the 
same way as explained for SU(2) and SU(3). 
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